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We present a general subtraction method for computing radiative corrections to QCD jet cross sections at 
next-to-next-to-leading order accuracy. The steps needed to set up this subtraction scheme are the same as 
those used in next-to-leading order computations. However, all steps need non-trivial modifications, which we 
implement such that that those can be defined at any order in perturbation theory. We give a status report 
of the implementation of the method to computing jet cross sections in electron-positron annihilation at the 
next-to-next-to-leading order accuracy. 



1. Introduction 

Accurate predictions of QCD jet cross sections 
require the computation of radiative corrections 
at least at next-to-leading order (NLO) accuracy, 
but in some cases also at higher order. The 
physical cases when computations at the next- 
to-next-to-leading order (NNLO) are important 
have been discussed extensively in the literature 
[1]. Although perturbation theory is expected to 
be a rather systematic procedure, this has not 
been reflected for many years in the computation 
of radiative corrections to QCD jet cross sections. 
The main reason for this is that the higher or- 
der corrections are sums of several contributions 
which are separately divergent in d = 4 space- 
time dimensions, only their sum is finite. Further- 
more, these contributions have different numbers 
of particles in the final state therefore, their com- 
bination is not straightforward. 



2. Subtraction methods at NLO accuracy 

The perturbative expansion of any jet cross sec- 



,LO 



,NLO 



tion can formally be written as a = a'""^ +a" 
^NNLO ^ Lg^ ^s consider e~^e^ -^ m jet pro- 

duction, when a^'-' is the integral of the fully ex- 
clusive Born cross section over the available phase 
space defined by the jet function J^, 



LO 



dcr„jJm 



d0„|A4(^)|V„ 



(1) 



The NLO correction is the sum of two contribu- 
tions. We have to consider the fully exclusive 
cross section dcr^ for producing m + 1 partons 
and the one-loop correction da^ to the produc- 
tion of m partons. 



,NLO 



dcr^J, 



m+1 



771 + 1 



dCT^J, 
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j dcp,^+,\M^^l^^Jm+l 



d</.„,2Rc(A^,W|A^(o))J„ 



(2) 



These two contributions are separately divergent 
in d = 4 dimensions although their sum is finite 
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for infrared safe observables. We assume that 
ultraviolet renormalization has been carried out, 
so the divergences are purely of infrared origin 
and are regularized by defining the integrals in 
d = 4 — 2e dimensions. 

There are several general methods of comput- 
ing the finite NLO correction. Most of these rely 
on the same principles, namely one defines an ap- 
proximate cross section da^ that regularizes the 
real correction in d dimensions in all its infrared 
singular limits, so the cross section 

^r+^ = / [{da^).-o Jrrr+i (da^),_„ J™] 

Jm+l 

(3) 

is finite. The subtraction term in this equation is 
symbolic in the sense that it is actually a sum of 
different terms and the jet function depends on 
different momenta in each of these terms. 

The first step in defining the approximate cross 
section is to derive the universal factorization 
properties of QCD matrix elements when one ex- 
ternal momentum becomes soft, or coUinear to 
another one (refered to as unresolved). These are 
well-known both at NLO and at NNLO. In writ- 
ing the factorization formulae we use the colour- 
state notation [2 and also some notation intro- 
duced in Ref. [3]o For the case when parton r 
becomes soft, we have 



Sr\Mt\M 



■)? 



i,k 



Hr^kr 



{M^°H---)\T^-Tk\M^^H---)) 



(4) 



and for partons i and r becoming collinear, we 
have 



C,r\M^^X,{p,,Pr,...)\'^ 




±{M^^Hp,^, . . .)\P^^\z„ kr, e)\Ml^{p.r, ■ 


■•)), 




(5) 



(0) 



where Sji ~ 2pj -pi (j, I — i, k, r), P^^' are the 
AltarcUi-Parisi splitting kernels and Zi denotes 



^We drop some numerical factors in order to keep the ex- 
pressions as simple as possible, as only the structure of 
these formulae is relevant for the discussion. 



the momentum fraction of parton i to the total 
momentum of the splitting parton. Similar for- 
mulae can be written for describing the IR struc- 
ture of the squared matrix element at NNLO. The 
necessary ingredients of these factorization for- 
mulae, namely (i) the tree level three-parton split- 
ting functions [4|5|6|7|8|9|10|lll and double soft 
gg and qq currents |ll|12j . (ii) the one- loop two- 
parton splitting functions |13I14I15I16| and soft- 
gluon current T7] have been known for some time. 
The difficulty of using the multiple infrared fac- 
torization formulae for constructing the approxi- 
mate cross sections is amply demonstrated by the 
slow progress in setting up a (general) subtraction 
scheme. 

The second step is to write the IR factorization 
formulae in such a way that intersecting limits 
can be identified and disentangled so that multi- 
ple subtraction is avoided. At the NLO accuracy, 
the only such intersection occurs in the regions of 
phase space where one parton is simultaneously 
soft and also collinear to a second (hard) parton 
and the overlap of the soft and collinear limits 
can easily be identified to be the collinear limit of 
the soft factorization formula [3j, 






(6) 



Sir J- ^i 

Thus the candidate subtraction 



E 









(7) 



has the same singular limit as the real correction 
itself and is free of double subtractions. How- 
ever, disentangling the unresolved limits at higher 
orders, when multiple soft, collinear and soft- 
coUinear limits overlap in a complicated way, is 
far more cumbersome 3^. This calls for a simple 
and systematic procedure. 

In a physical gauge the collinear singularities 
are due to the collinear splitting of an external 
parton [18119] . The overall colour structure of 
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the event does not change, the sphtting is en- 
tirely described by the AltareUi-Parisi functions 
which are products of colour factors and kine- 
matical functions describing the dynamics of the 
coUinear splitting. The emission of soft gluons 
is just the opposite; it does not affect the kine- 
matics of the radiating partons, but it does affect 
their colour, because it always carries away some 
colour charge. If we want to identify the coUinear 
contributions in the soft factorization formulae to 
any order in perturbation theory, we can use the 
following simple procedure: (i) employ the soft 
insertion rules [20|llj to obtain the usual expres- 
sion 



1=1 fc=l hoi. 



with 



fiiy 



pi^n" 



■p';:n'' 



hcl. 



Pr-n 



(8) 



; (9) 



(ii) fix the gauge vector to n^ — Q^ 



nMn2 



/SrQ, 



SrQ = 2pr-Q {Q^ is the total incoming momen- 
tum) to identify the coUinear contribution in the 
colour-diagonal terms 






E 



2SfcQ 

._.. ■^ ... \ SirSrk SrQSir SrQSkr 



JE 



x{M^^K■■■)\T^■Tk\M'^:^K...)) 



(0), 






(10) 



(iii) define momentum fractions in the Sudakov 
parametrization of momenta p^ and p{f being 

-, so that the colour- 



collinear as z,- = — "- , 
diagonal terms become equal to the coUinear limit 
of the soft factorization formula. Then the pure 
soft contributions are given by 



E 



- m 



2s. 



ik 



ZSj, 



2s 



kQ 



L \ SirSrk ^rQ^ir ^rQ^kr 



(0)/ 



x{Mlli>i...)\T,-T,\Mi::> {...)) 



(11) 

We checked explicitly that this procedure leads to 
non-overlapping factorization formulae that de- 
scribe the analytic behaviour of the squared ma- 
trix elements in any IR limit at the NNLO accu- 
racy [21] . 

The third step is the definition of the subtrac- 
tion terms. The factorization formulae in Eqs. ^ 
and (O (and similar ones at NNLO) are valid in 
the strict limits only and have to be extended over 
the whole phase space. This extension requires 
momentum mappings {p}m+i ~^ {p}m that 

• implement exact momentum conservation, 

• lead to exact phase-space factorization, 

• can be generalized to any number of unre- 
solved partons and 

• respect the (delicate) structure of cancella- 
tions among the subtraction terms. 

For any such mapping we can write the exact fac- 
torization of the phase space in the following sym- 
bolic form 



d0™+i(Pi 



;g) = d0„,(pi,...;Q)[dp,]. (12) 



Then the singular integral over the momentum of 
the unresolved parton r can be computed inde- 
pendently of the jet function and the rest of the 
phase-space integration, leading to 



f da^ = da^ J(e) , 



(13) 



where /(e) is an operator in colour space with 
universal pole part. 



i 

+0(e") , 



^'-^£--(^)' 



(14) 
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with 7i being the flavour constants defined in for 
instance, Ref. P|. 

The fourth step is to identify the universal IR 
pole structure of one-loop QCD matrix elements, 

\M^1H{P})) = -ll{e)\Mi^H{p})) + Oie"). (15) 

We observe that the poles in the insertion oper- 
ator 7(e) are equal, but opposite in sign to the 
poles of the virtual correction, so that the m- 
parton integral 



,NLO 



da' 



da' 



(16) 



6 = 



is finite and we can take the physical limit e ^ 0. 
Therefore, the sum of the two finite contributions 
a^'^ and a^^ can be computed numerically in 
four dimensions and is equal to a^^'^ . 

3. A subtraction scheme at NNLO accu- 
racy 

The physical motivation for higher accuracy 
and the success of the subtraction schemes at 
NLO lead one to consider the extension of the 
subtraction method to NNLO, when three terms 
contribute: the double-real, the real-virtual and 
the double- virtual cross sections, 



,NNLO 



RR 

'm+2 



RV 

'm+1 



VV 



/ dcr,„_^2'^m+2 + / dcr^+i J,„+i , 

Jm+2 Jm+1 y 



7) 



J VV T 



The reorganization of the NNLO contributions 
into three finite cross sections, 



^NNLO 



NNLO 
'm+2 



NNLO 
•'m+1 



^NNLO 



is governed by the jet function as follows: 



^NNLO 
'^m+2 



^NNLO 
'm+1 



m+2 
m+2 



{da, 



RR T 

+2-Jm+2 



_ , RR.Aa J 



1 RR,A^ J 
C10'^_Lo Jm+1 



icr, 



^ {( 

m+r- ^ 



da. 



RV 

m+1 



do-, 



RV,Ai 

m+1 



da. 



_RR,A 

m+2 

J RR,A, 
dCT™+2 

RR,A, 

m+2 



•Jm j f J 



J, 



rn+1 



j •Jm f J 



(18) 



(19) 



(20) 



.NNLO 



{d-r+/(dC2^^-da, 



RR,Ai: 
m+2 



da. 



RV,A, 



^ + {l^-T+tY^]}jrr 



(21) 



Here we see that at NNLO accuracy one has 
to disentangle the overlapping singularities also 
among the singly- and doubly-unresolved limits. 
The purpose of the approximate cross section 



do'm+2 ^^ i^ to regularize the singly-unresolved 

RR A 

limits of dCT^^'2 ^ and the doubly- unresolved lim- 
its of dcr„j^2 ^ simultaneously. Similarly, the ap- 
proximate cross section I ^^ dCT„j^'2 ^ \^^ regular- 
izes the singly- unresolved limits of j^ '^<^m+2 ^ ^'^d 
the e poles of dcr„j^\ ^ , respectively. This puts 
severe constraints on the phase-space mappings 
needed for the definition of the subtractions. 

In Ref. [231 we introduced new types of momen- 
tum mappings, one for coUinear- and another for 
soft-type subtractions, that can easily be general- 
ized to any order in perturbation theory. The key 
feature of these mappings is that in the factorized 
TO-particle phase space all momenta take away the 
recoil instead of a single one as in the case of the 
dipole (or antennae) subtractions. In this way the 
factorization of the phase space can be done in a 
way that respects the (delicate) structure of can- 
cellations among the various subtraction terms. 

The complete subtraction scheme at NNLO, 
based on these new, fully local approximate cross 
sections is defined in Refs. |22|^3124| . We em- 
ployed this subtraction scheme for computing the 
finite cross sections a^^i^ and af^^i^ of the C- 
parameter and thrust distributions in electron- 
positron annihilation. Fig.[T] shows the distri- 
butions normalised to the total cross section at 
0(as) accuracy. The computer time needed for 
obtaining these distributions is fairly little. The 
plots shown here can be obtained on a desktop 
computer in about 50 hours. The still missing 
three-parton contribution is a smooth function as 
compared to the four- and five-parton contribu- 
tions, therefore its numerical integration does not 
raise any serious stability issues. 
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Figure 1. The five- and four-parton contributions (lower panels) to the C-parameter (left) and thrust 
(right) distributions. The upper panels show the distributions at the LO, NLO and incomplete (without 
the three-parton contribution) NNLO accuracy. 



In order to have the complete physical predic- 
tion we also have to add a^^'~', which requires 
the integration of the subtraction terms over the 
singly- and doubly-unresolved factorized phase 
spaces. 

The necessary one-particle integrals have been 
computed in Refs. |25l26l27j . In Ref. g5] we used 
standard techniques of partial fractioning, iter- 
ated sector decomposition |28|29j and residuum 
subtraction to find the Laurent expansion of the 
one-particle integrals in 



da: 



RV,Ai 

1 



da 



RR.Aj^Ai 
m+2 



(22) 



With this technique the expansion coefficients are 
given as finite integrals. In order to find explicit 
analytic expressions, one can use the integration- 
by-parts technique to reduce the integrations to 
master integrals to be evaluated by solving differ- 
ential equations pB' . An alternative solution is to 
employ the Mellin-Barnes technique to evaluate 
the expansion coefficients in terms of harmonic 



sums [57]. Our study shows that the latter tech- 
nique is in general more efficient. We expect that 
the same techniques can also be employed for the 
computation of the coefficients in the e-expansion 
of the two-particle integral 



da. 



RR.Ag 

m+2 ' 



da. 



RR^Ai; 



(23) 



This work is in progress. 

In this contribution we have discussed how the 
same steps of setting up a subtraction scheme for 
computing QCD jet cross sections at NLO accu- 
racy can be used, but need to be modified when 
generalizing to NNLO accuracy. We introduced 
these modifications such that a systematic pertur- 
bative expansion can in principle be given at any 
order of perturbation theory. The explicit com- 
putations of course become rather cumbersome 
already at NNLO. 

We are grateful to our collaborators: U. Agli- 
etti, P. Bolzoni, V. Del Duca, C. Duhr and 
S. Moch. 
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